o 

(N 



Oh- 






The Gough-James Theory of Quantum Feedback 
Networks in the Belavkin Representation 

O. G. Smolyanov 1 , A. Truman 2 



Abstract 

The mathematical theory of quantum feedback networks has recently 
(f) ' been developed by Gough and James [5j for general open quantum dy- 

namical systems interacting with bosonic input fields. In this article we 
show, that their feedback reduction formula for the coefficients of the 
closed-loop quantum stochastic differential equation can be formulated in 
terms of Belavkin matrices. We show that the reduction formula leads to 
-—t a non-commutative Mobius transformation based on Belavkin matrices, 

and establish a *-unitary version of the Siegel identities. 
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1 Introduction 

l/~) ■ In a recent publication [5] Gough and James have introduced a model for a 

quantum feedback network. Each component may be modelled in isolation as 

f***^ ■ a Hudson-Parthasarathy model, valid for quantum optical models, and repre- 

sented as a single vertex with an equal number of external inputs and outputs 
carried along semi- infinite transmission lines represented as directed edges. The 
algebraic procedure is to collect all the operator coefficients of the associated 
quantum stochastic differential equation governing all components into a ma- 
trix. This gives the open-loop description, and feedback is introduced by con- 
C^ ■ necting various input/output lines to giver internal edges. They use a Hamilto- 

nian model for the entire network which generalizes the Chebotarev-Gregoratti 
Hamiltonian which describes the propagation of the fields along the edges and 
their interaction at the vertices. They obtain a Markovian limit for the network 
in a zero time delay limit, eliminating the internal edges in the process. The 
limit quantum stochastic differential equation is then described by a reduced 
operator matix. 

We show that their feedback reduction formula is a Mobius transformation 
that can naturally be extended to mappings into the Belavkin matrix represen- 
tation of quantum stochastic calculus. In particular we interpret this as a non- 
commutative Mobius transformation between *-unitaries. We begin by recalling 
the basic notions of quantum stochastic calculus and its Belavkin formulation. 



1.1 Quantum Stochastic Processes 

The Hudson-Parthasarathy theory of quantum stochastic calculus considers 
quantum stochastic processes as operator valued processes on Hilbcrt spaces 
of the form Sj = f) ® T (J <g) L 2 [0, oo)) where F)o is a fixed Hilbert space, called 
the initial space, and t is a fixed Hilbcrt space called the internal space. We 
shall be interested in the finite dimensional case 6 = C™ where n > 1. Here 
r (t) denotes the second quantization functor to (Bosonic) Fock space. We shall 
denote the time variable t by A® . Taking {\ei) : i = 1, • • ■ ,n} to be an or- 
thonormal basis for t, the creation process to state \ei) will be denoted as Af, 
while its adjoint, the annihilator for the state, is denoted as A® 1 . The scattering 
process from state \ej) to state \ei) will be denoted as A] 1 . In this way, we 
have the (1 + n) x (1 + n) fundamental quantum processes Af . (We adopt the 
convention that Latin indices range over 1, • • • , n while Greek indices range over 
0, 1, • • • , n. We also apply a summation convention for repeated indices over the 

corresponding ranges.) We note that we have I A" ) = A^ . 

As is well known S) decomposes as -f)[o,t] <8> $)(t,oo) f° r each t > where 
Sj [ot] = f) (gi r(t®L 2 [0,£)) and % :CO ) = r (t <g> L 2 (t, oo)). We shall write 
2l t i for the space of operators on S) that act trivially on the future component 
fi(t,oo)- A quantum stochastic process X t = {X t : t > 0} is said to be adapted 
if X t <G 2l t] for each t > 0. 

Taking {x Q ^ (t) : £ > 0} to be a family of adapted quantum stochastic pro- 
cesses, we may then form their quantum stochastic integral X t — L x a p (s) dA" 13 
where the differentials are understood in the Ito sense. Given a similar quantum 
Itb integral YJ , with dYt — y a /3 (t) dA^ , we have the quantum Ito product rule 

d (X t .Y t ) - dX t .Y t + X t .dY t + dX t .dY u (1) 

with the Ito correction given by 

dX t .dY t = x ak (t)y k[i (t)dAf. (2) 

The coefficients {x a p (£)} may be assembled into a matrix 



x = j a;oo (t) 



X,Q (t) 



x 0» (*) \ e a (l+n)x(l+n) /gx 



X.. (t) J c ~*1 



which we call the ltd matrix for the process. (Here we use the convention 
that xo,(t) denotes the row vector with entries (xoj (t))™ =1 , etc. The Ito 
matrix for a product XtYt of quantum Ito integrals will then have entries 
{x a p (t) Y t + X t y aj3 (t) + x ak (t) y k[j (<)} and is therefore given by X t Y t +X t Y t + 



X t PY t , where P := 







1.2 Belavkin's Matrix Representation 

We consider the mapping from Ito matrices X £ gj(i+n.)x(i+n) ^ Q assoc i a ted 
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where I„ is the nxn identity matrix. The subscripts n will be generally dropped 
from now on for convenience. We have the following identifications 



X T 

XPY 

XY 



X* := JX f J[, 

XY, 

XJY. 



We shall refer to X* as the *— involution of '. 
x a/ 3 (t) dAf may then be written as 



The Ito differential dX t 



dX t =tr{x t dA t \, 
where (with / denoting the usual transpose for arrays) 







dA t := (dA ') 



dA 00 







{dA")' 
(dA'°)' 
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The main advantage of using this representation is that the Ito correction XPY 
ca now be given as just the ordinary product XY of the Belavkin matrices. 

Let X t and Y t be quantum stochastic integrals, then the quantum Ito product 
rule may be written as 

d (X t Y t ) =tr{ [(X t I + X t ) (Y t I + Y t ) - (X t Y t ) I] dk t } . (5) 

The process / (X t ) has differential df (X t ) = tr{[f (X t I + X t ) - f (X t ) I]dA t }. 

1.3 Evolutions and Dynamical Flows 

Hudson and Parthasarathy [11] show that the quantum stochastic differential 
equation (QSDE) 

dV t = tr [GVtdAtj , Vb = l, (6) 



has a unique solution for a given constant Belavkin matrix G = V — I of co- 
efficients on 03 (f)o), the bounded operators on f) . Necessary and sufficient 
conditions for unitarity are then given in Belavkin representation by 

(I + G) (I + 6)* = I = (1 + G)* (I + G) . 

This states that V = I + G is ^-unitary on 03 (t)o) (1+ " +1)x (1+n+1) , that is 

W* = I = V*V. (7) 

We may write the QSDE as dV t = tr j [YV t - W t ] dA t \. 
Lemma 1 The most general form for V leading to a unitary is 
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(8) 



^O'O' / 



where S is a unitary in 03 (f)o) , L is a column vector length n with entries 
in 03 (f)o)> an d H is self-adjoint in 03 (f)o)- 

The proof follows from the analysis of [TT]. The triple (S, L, H) is termed the 
Hudson-Parthasarathy parameters of the open system evolution. In standard 
notation the QSDE reads as (summ over all field multiplicities) 



dV t 



{{S jk - 6 jk ) dAl k + L 3 dA? - L*S jk dA° t j - (h*^ 



iH dAS 



V t . 



We interpret A^ as the input noise and V t *A"Vt as the output noise. 



2 Quantum Cascaded Systems 

If two systems are cascaded in series then the Hudson-Parthasarathy parameters 
of the composite system were shown to be [B],[5] 



o St 
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Here the output of the first sytem (Si,L\,Hi) is fed forward as the input to 
the second system (52, £2, H 2 ) and the limit of zero time delay is assumed. As 
remarked in [5], the series product actually arises natural in Belavkin matrix 
form as 



V 



V 2 Vi. 



The product is clearly associative, as one would expect physically, and the 
general rule for several systems in series is then V ser i e s = V n ■ • • V2V1 



3 General Feedback Reduction Formula 



The internal edges may be eliminated in a zero time delay limit to obtain a 
reduced model. Let < n± < n be then number of internal edges to be elimi- 
nated, and let n e = n — n± be the remaining edges. The algebraic information 
about the original network is contained in the matrix V which we partition as 



V = 



Here we decompose indices into two groups e and i distinguishing external and 
internal. That is, Voo = Vq>o' = l,V e0 = V i0 = V >o = V > e = Vo'i = 
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and [V 0e V 0i ] = -SL*. 

Theorem 2 We assemble a Belavkin matrix T(V,X) in 2l(i+«.+i)x(i+™.+i) 
with sub-blocks 

T (V, X) a(j = V al3 + V al X (1 - V^xy 1 V ip 

for a — 0, e, 0' and (3 = 0, e, 0', where we fix a unitary operator in X £ C" lXni 
such that the inverse above exists. Then J- (V, X) is again a ^-unitary, that is, 



T (V, X)* T (V, X) = T (V, X) T (V, X)* = I, 



(9) 



so that T (V, X) determines a unitary dynamics for the reduced set of n e inputs. 
Moreover, we have the identity 



F(V,X)* = F(V*,Xi) 



(10) 



Remark 3 The matrix X appearring above is typically an adjacency matrix 
in applications, describing which internal outputs are to be connected to which 
internal inputs. In engineering, it could be interpreted as a gain matrix. We 
also point out that the involutions in (|j^p are on spaces of different dimensions. 
The first involves J ni while the second involves J„. 

Proof. The construction of T (V, X) is essentially the rephrasing of the 
Mobius transformation associated with the reduction, introduced in [5], in the 
language of Belavkin matrices. The construction in flS) clearly yields a Belavkin 
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matrix over the remaining n e external degrees of freedom. It is a straightforward 
calculation to show that, for X unitary, the matrix takes the form 



F(V,X) 



with the Hudson-Parthasarathy parametrizing operators (iS , L rcd , _ff rcd ) given 
by 

L = L e 4- S ei (X — Sa) Li, 

# rcd = tf+^ImLJS-^X- 1 -^)^, 

2— i,e 

in agreement with [5]. ■ 

For X e 2l( 1 +"+ 1 )x( 1 + m + 1 ) we shall introduce the extended convention X* := 
JJ m XTJ„. Let V be the Belavkin matrix generating a unitary quantum dynamics 
as above, we define the Mobius transformation $ : D t-> 2(( 1 + n .+ 1 )x( 1 + n «+ 1 ) by 
$ = J"(V, •) with domain D = {X e C™ 1 *™ 1 : I - V ti X is invertible}. 

Theorem 4 TTie mapping $ satisfies the Siegel type identities 

$(xf$(y) = i+ y el (l-x^l) (xty - i) (i - ^y)- 1 

$(X)$(F)* = l + {V i0 ,V lj .,V i sy)(l-XV 1 . 1 )- 1 (XY^-l)(l-V^Y) (W,, r le . r i() ,)* . 
In particular, $ maps unitaries to -k-unitaries. 

Proof. The form of these relations are similar to the standard Siegel iden- 
tities, see for instance |14) . but with the ^-involution now replacing the usual f. 
The algebraic manipulations involved are otherwise identical. ■ 

We remark that the standard Siegel type identities have independently been 
extended in an entirely different direction to deal with Bogoliubov transforma- 
tions in a recent paper of Gough, James and Nurdin [15J . They replace the usual 
f-involution with an alternative involution, this time on the space of doubled up 
matrices required to describe the symplectic structure, however they similarly 
rely on the argument used in the above proof. 
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